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ABSTRACT 


The Heat and Mass Transport of Furniture Fire, Smoke and 
Toxic Gases (HEMFAST) program simulated the furniture fire and the 
Lomc fas, growth un’ a room and its spread to the other rooms. § This 
was achieved by incorporating an improved furniture fire model into 
uge FAST model in current. usage iat NBS. The furniture fire model 
was designed to be self-controllable in Simulating a fire source 
utilizing as inputs only the bench scale data of new scalable 
material burning rates, flame spreading, thermal ignition, and the 
flame properties of the following: soot volume fraction, C50 and 
H50 mass fractions of fuel weight, and gas temperatures. Since the 
Hottels' zoning method was used in the thermal radiation model it 
was possible to consider surface emissivities less than unity on 
all surface elements, and to consider three isothermal gray gas 
volumes which are the high-temperature flaming plume, the sooty 
upper gas layer, and the lower gas layer at 50% relative humidity. 
Some salient features are noted in the simulation of a sample fire 


growth on a furniture mockup. 
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SECTION 1 
INTRODUCTION 


A milestone has been reached in the development of a model 
which simulates the furniture fire and the toxic gas growth ina 
room and its spread to other rooms. This was achieved by merging 
an improved version of the Furniture Fire Model with the FAST 
program, a current multi-room fire model developed at NBS. The 
result is the Heat and Mass Transport of Furniture Fire, Smoke, and 
Toxic Gases (HEMFAST) computer program. The Furniture Fire Model 
has three important qualitative programming features. The first is 
that each main module was designed to be flexible, efficient, and 
modular. Secondly, each main module has made effective use of bench 
scale data such as the ignition and flame spread rate, the burning 
and constituent mass loss rates, and the thermal radiation 
properties of furniture samples. Third, each main module was able 
to fit into the specific areas in the FAST code, making it more 
likely that the HEMFAST modules will fit in other room fire models, 
such as the HARVARD CFC series. A high level flowchart for HEMFAST 
is shown in Appendix A and is an elaboration on these comments for 
specific modules. 


The main modules of the Furniture Fire Model were described 
in the interim report, "Furniture Fire Model", Reference 1, and 
their high level flowcharts are reproduced in Appendix A. The 
Furniture Fire Model Demonstrator was described in the next interim 
report, "Mathematical Modeling of Furniture Fires", Reference 2, 
and was improved to obtain more realistic predictions of the fire 
growth and to prepare the model for merging with FAST. In the 
interim report a new model of thermal ignition and flame spreading 
was described and was applied to the demonstrator. In this work 
the application of the thermal ignition and the flame spreading 
models was refined to more accurately calculate the pre-ignition 
parameters such as the preheat temperature and the pre-ignition 
heat flux distribution located from the flame front out to the 


effective location of the preheat temperature. The module for the 


Local Burn Rate, Fuel Consumption, and Convective Heat Flux 
Predictors (see Appendix A) was improved by considering the effect 
of net. total surface heat fluxes on the local burn rdate andes 
fuel consumption. This required the development of new scaling 
relationships for the data from the cone calorimeter for applica- 
tion to the furniture fire modeling. In the Flame Structure Model 
module, the conical frustrum was replaced by a polygonal based 
parallepiped to more accurately represent the flame shape observed 
in video recordings of furniture fires conducted in the NBS 
Laboratory (Reference 3). The flame height modeling was also 
Slightly modified. The approach to calculating the flame radiative 
properties is that of Bard and Pagni (Reference 4). 


Examination of Version 17 of the FAST code and discussion 
with the NBS Center of Fire Research personnel indicated the FAST 
code needed two major improvements. The first was to insert a 
self-consistent fire source model in which the fire growth responds 
dynamically to its environment: The°main response “factors Char 
need to be taken into account are the radiative feedback for con- 
trolling burn rate and flame spread, the air entrainment into the 
plume, and the air constituents concentration. The second major 
improvement was the extension of the Thermal Radiation Heat 
Transfer Analyzer module in the furniture fire model to a more 
detailed radiative heat transfer analysis of the room, including 
the effects of soot, carbon dioxide, and water vapor in the upper 
gas layer and of the water vapor in the lower gas layer. The air 
entrainment modeling seemed adequate in the FAST code itself. The 
Only response factor not yet considered was that of reduced oxygen 
concentration in the lower gas layer as a result of its mixing with 
the upper gas layer. This response factor awaits the availability 
of the bench scale data of flame spread and burn rate due to vary- 
ing ambient oxygen concentrations. 


The resulting computer ‘code is HEMFAST. Its*> input* and 
results are described in the’ last section ofthis report?*" In 
particular, a generic situation of fire growth on a piece of furni- 


ture in a room with one opened window was reasonably simulated. 


The furniture burnrate, the temperature profiles of the room and 
the gas layers, and the upper gas layer constituents concentrations 
are displayed as a function of time in Section 3.3 of this report. 


These predictions are qualitatively consistent with the furniture 
mockup tests in Reference 3. 


SECTION 2 
FURNITURE FIRE MODEL IMPROVEMENTS 


The purpose of the improvements was to add more realism to 
the demonstrator version and to prepare the furniture fire model 
for incorporation to the FAST program. There are four main im- 
provements covered in this section, which are: (1) the fine-tuning 
of the flame spread algorithms, (2) new computations of the local 
burn rate in response to varying surface heat flux, (3) a new flame 
shape algorithm, and (4) incorporation of radiative properties by 
Bard and Pagni. The modification of the thermal radiation analyzer 
module to take into account the room and the gas layers geometries 
is explaineda lastly. 


Zia FLAME SPREADING ALGORITHM IMPROVEMENTS 


The flame spreading algorithm was improved by incorporating 
better interpretations of the flame spread rate equations in 
Reference 2. The flame spread rate formulae for thermally thick 
and thin materials are reproduced in a slightly modified form 
below. 


Thermad lyathnickvyd imat. PO 4% * povyuyso; for RMMA: 


2 
2B 
ble v4 S/T 
¥vpck ie. 1 


Thermally thin limit, peck +0, for polyurethanes: 
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That is, the preheat temperature, Tas is no longer the surface 
temperature at ignition time due to the net surface heat fluxes 
other than the flame feedback heat flux. i is now defined as the 


Current surface temperature at some surface distance, So» from the 
flame front on the virgin material. This implied a different 
interpretation of the feedback heat fluxes to mean the net surface 
heat flux distribution between the flame front and the surface 
distance, Sa: Dns near tiluxjdistribucionsis used insveither flame 
Spread rate formulae listed above by constructing qr, and S, as 
shown in Figure 1, and using the appropriate summation. If the 
heat flux distribution is some well defined function of s, the 
Summation can be converted to a convenient integral. This new 
interpretation of the flame spread rate formulae is significantly 
different than commonly encountered in the literature, but it was 
necessitated by the requirement for ease of application to inodeling 
flame spreading in the furniture fire model. 


The flame spread procedure used in the furniture fire model 
foetuses vraved min Figure 2. A furniture cushion jpariel is divided 
up into square-like surface elements and nine of elements are shown 
in the figure. The net surface heat flux and the surface tempera- 
ture is computed as a function of time at the midpoint of each 
element. Thermal ignition is assumed to occur when the surface 
element; vemperature reaches the ignition temperature. Piloted 
ignition is assumed to occur when the program user selects the time 
maiomechne; Location of ignition. -Once ignition occurs at a location, 
a small octagonal burning area is constructed. Flame spreading is 
simulated by moving the vertices of the burning polygon in a normal 
outward direction over an incremental distance in an incremental 
mime per1o0.- Ine preheat tenperature, Tis the preheat heat flux, 
i? corresponding to the 
burning polygon vertice, i, are selected in the following way. 


qi» and the preheat surface distance, s 


-Each virgin surface element temperature and heat flux is extrapo- 
lated to the intersection of the flame spread direction vector and 
of the vector drawn through the surface element midpoint parallel 
to the appropriate polygon line segment, as illustrated in Figure 
2. Several such intersections can occur on a single flame spread 
vector because there are several virgin surface elements. The 
intersection point closest to the vertice of the flame spread 


Figure l. Construction of q, and s; from Heat Flux 
Distribution, q = f(s), Fixed Ahead of the 
Flame Front at Speed Vp° 


Burning 


+ 


Polygon 


Figure 2: HEMFAST Flame Spread Procedure on a Heated 
Cushion Panel. 


direction vector then determines the preheat surface distance, the 
preheat temperature, and the preheat heat flux....The.closest burn- 
ing surface element|to the flame*spread\direction veclormaa, | 
Similarly be found and the net surface heat flux of the chosen 
burning surface element is extrapolated to the burning polygon 
VELCUCEX. 


Note the net surface heat fluxes can be decomposed into the 
radiative and the convective portions. The pre-ignitive heat flux 
distribution .required forythe flame spread. rate formulac saeco. 
structed by fitting an exponential function between the radiative 
heat flux at the polygon vertex and the preheat radiative heat 
flux, and fitting a step;function for the convective Neate: 
distribution. The convective heat flux distribution fseedgtTevem. 
for opposed or assisting airflow, and so the appropriate boundary 
layer theory must be used. Other heat flux terms not accounted 
for, such as the flame foot convective/conductive heat flux or the 
glowing combustion/oxidative pyrolysis feedback heat flux, are 
treated as a constant to be curve fitted to the flame spread bench 
scale data. As a last comment, the flame spreading algorithm 
includes features to merge overlapping burning polygons, to keep 
Spacings between polygon vertices at reasonable values, and to 
initiate flame spreading into adjoining cushion panels. An example 
of flame spreading simulated on a generic chair with a four cushion 
panel configuration divided into a total of 64 surface elements is 
shown in Figures 3 to 6. 


ead NEW SCALING OF THE CONE CALORIMETER DATA 


The Furniture Fire Model Demonstrator does not adjust the 
burning rate data from the cone calorimeter to account for varying 
irradiance levels, because the particular data used was only col- 
lected ata Single irradiance. lever of 2. kW/m. Also, there was 
no effective way to scale the burning rate data as a result of the 
varying net surface heat flux. In an effort to develop such a 
Scaling relationship for the furniture materials tested in the cone 


calorimeter, a large series of the cone calorimeter database in the 
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Example of Flame Spreading on Mockup Seat at 20, 40, 


80, and 100 Seconds on Surface Elements 1 to 16. 
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RIGHT ARM WIDTH (M) 


Example of Flame Spreading on Right Armrest at 40, 60, 

80, and 100 Seconds Simulation Time on Surface Elements 17 
to 32. (This armrest begins burning near the center at 40 
seconds and the flame spreads toward both edges. The 
flame spreading on the armrest is driven by the growth of 
the burning polygons on the seat cushion as these polygons 
intersect the armrest cushion. The dark lines in the 
figure show the flame positions at the times provided.) 
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BACK WIDTH (M) 


Example of Flame Spreading on Backrest at 40, 60, 80, and 
100 Seconds Simulation Time on Surface Elemens 33 to 48. 
(This backrest begins burning near the center at 40 
seconds and the flame spreads towards both edges. Flame 
spreading on the backrest is driven by the growth of the 
burning polygons on the seat cushion as these polygons 
intersect the backrest. cushion. The dark lines in the 
figure show the flame positions at the times provided.) 
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LEFT ARM WIDTH (M) 


Example of Flame Spreading on Left Armrest at 80 and 100 
Seconds Simulation Time on Surface Elements 49 to 64. 
(This armrest begins burning near the left side at 80 
seconds and the flame spreads toward both edges. Flame 
Spreading on the armrest is driven by the growth of burn- 
ing polygons on the seat cushion as these polygons 
intersect the armrest cushion. The dark lines in the 
Figure show the flame positions at 80 and 100 seconds.) 
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form of tabular output were obtained from John Krasny and Vito 
Babrauska (see References 3 and 5). This series included data 
corresponding to the following furniture materials: cotton 
fabric/treated PU foam, cotton fabric/non-treated PU foam, 
polyolefin fabric/treated PU foam, and polyolefin fabric/non- 
treated PU foam. These materials were subjected to cone 
calorimeter irradiances of 25, 30, 40, and 50 kW/m* during the 
tests. These 16 experimental runs were used to develop scaling 
relationships that allow us to use the concepts of heat of 


pyrolysis and the heat of combustion in the furniture fire model. 


The scaling relationships obtained are shown in the following. 


Sealed heat release rate and mass loss rate: 


.* hs i 
Deis : (3) 
qj Vo x 
o* Qo e 
Oh epg ly (4) 
(=. | a 
where 
4 
Ae eanQawt qr Rio, 9 Ty + q, = net surface heat flux (5) 


Sealed time is best defined by the heat release: 


ae es) nj Celie G,) dt, 
Oo 
where 
* i ° 
Q - Q; es nh; dt 
G7. = cos | 2m Smh ie + 1 
i ay Q 
i { 
4 n 
iene 5 Q,/n = averaged total heat release of the (8) 


test sample i 


ys 


corresponding with h,and m,- 


(7) 


(9) 


The furniture materials are known to undergo glowing combus- 
tion or smoldering at low external irradiances, even when there are 
net surface heat flux losses. Here the net surface heat flux, Qj» 
is defined as the sum of. the surface convective,. the flame) surface 
irradiant, the material surface radiant, and the external irradiant 
heat fluxes. These four surface heat fluxes can readily be calcu- 
lated and their formulae are the same as the ones, used in the 
furniture fire model. The surface convective heat flux, do: is 
calculated in the pool fire using the correlations by Orloff and 
deRis (Reference 6) and in the wall fire using the integral bound- 
ary layer model by Admah and Faeth (Reference 7). In both cases 
the surface convective heat flux is primarily sensitive to the 
material local mass loss rate, to the flame temperature and to the 
surface orientation. The flame surface irradiant heat flux, dnp? 
is calculated using the mean beam length theory but with the same 
flame structure and the flame radiative properties as used in the 
furniture fire model. It is noted this irradiant Neav {iui 
primarily sensitive to the burning rate, flame geometry and tem- 
perature, and to the maximum soot volume fraction in the flame. 
Since the burning material can be treated as a black-body (to a 
good approximation), the material radiant heat flux is assumed 
equal to the Boltzmann's constant times the fourth power of the 
burning temperature. The external irradiant flux, dri? i832 
measurable quantity in the cone calorimeter corresponding to the 
four irradiance levels. The surface heat flux, do? from the glow- 
ing combustion at low irradiance levels (or accounting for an 
inaccurate burning temperature) is very difficult to measure or 
calculate, and so is treated as a fitting constant and. isnusedmace 
basis to scale the heat release rate, he and fuel loss rate, m, to 
the starred ordinates, foe and, ae respectively as indicated in 
above equations. 
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A term for the material surface conductive heat flux is not 
included because one of the purposes of scaling the burnrate data 
is to avoid a heat and mass transfer analysis within the material. 
For a surface burning at a constant rate on an ideal thermally 
thick material, a theoretical heat conduction solution with phase 
change gives the result that an imposed surface heat flux divided 
by the fuel loss flux is equal to the sum of the heat of pyrolysis 
and the material enthalpy increase up to the burning temperature. 
Whereas for a thermally thin material whose temperature is already 
near the burning temperature, the material enthalpy increase need 
not be considered along with the heat of pyrolysis. Thus fora 
thermally thick PMMA material and for a thermally thin paper the 
concepts of heat of pyrolysis and the heat of combustion are useful 
meesocatine co any Size fire. This is the motdivation for creating 
a scaled heat release rate, nee and the scaled fuel loss rate, ma, 
in the hope of collapsing the data onto a single point or a single 
curve. 


For a composite furniture material that also can undergo 
charring, the effective heat of pyrolysis and the heat of combus- 
tion will change with the history of burning. Thus the scaled heat 
release rate and the scaled fuel loss rate will also change with 
the history of burning. Several researchers who have examined the 
curves of the heat release rate versus time for the different 
irradiances have attempted with some success to rescale the data as 
the heat release rate versus another abscissa variable, such as the 
amount of heat absorbed in the material, or the fuel loss, or the 


Nneaterebease asva funetion-of time. 


The variable corresponding to the heat absorption in the 
material was judged not applicable to this problem because of the 
additional complication of including the fitting parameter, do: 

The fuel loss variable was also avoided because it seemed too 
sensitive to extra bulk materials not involved in combustion such 
as moisture and other vaporizing inert material. The heat release 
as a function of time is a more tolerant variable to the con- 
taminants in the material. Also the total heat release of the same 


mo 


material under differing conditions of irradiances and humidity 
tend more to be of a constant than the total heat absorbed or the 
total fuel loss. The furniture materials however cannot be con- 
structed identically to each other, and furthermore, in real 
applications their foam thickness varies. Thus an additional 
scaling transformation must be applied to the heat release variable 
in order for the start-burn “and the “burnout *points to be latirne 
same value for similarly constructed materials under exposure to 
varying irradiances. 


One approach would be to normalize the time changing heat 
release of “the “maverial ie DY dividuals. wc oy Qi the total heat 
release at burnout time and multiplying by Q , the averagege ora. 
heat release of similar materials. This effectively stretches (or 
Shrinks) uniformly the time increments in the integral of the heat 
release rate. From intuitive considerations, however, there should 
be no stretching (or shrinking) of the incremental time intervals 
during the beginning phase of burning because the material behaves 
thermally like an infinitely thick translucent material with phase 
change, and also during the ending phase of burning because by then 
the material fs + thermally+thin. \ ‘Thus’ most of ‘thesstretchingsas 
time intervals should occur about halfway through the burning 
history. This is the reason for the cosine function of the normal- 
ized*héat ‘release "history in®the expression®for G. in the above 
equations. The amplitude of G; is defined such bid the scaled 
eet release (or scaled time or burn history), Tis is the same as 
Q at burnout ‘time for a constant value of the heat’ release ‘rate. 
Since the heat release rate is not a constant, the scaled heat 
release at. burnout..time,is.actually,a close approximations toe ne 
averaged total heat release of similar materials. 


The plots of the scaled heat release rate versus the scaled 
time and of the scaled fuel loss rate versus scaled time are shown 
in Figures 7 to 14 for the four different furniture materials at 
four different irradiances. The optimum value for the fitting 
constant, do? was 20 kw/m* for all four materials. All materials 


had in common the polyurethane foam, either treated or non-treated, 
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and either with cotton fabric or polyolefin fabric. The fact the 
foam can still sustain a fuel loss under glowing combustion condi- 
tions even with a net surface heat flux loss, gives credence to the 
idea that do is mainly the surface heat flux due to the glowing 
combustion or oxidative pyrolysis. We thus have a feature, not 
available in several published scaling relationships, for a smooth 
transition from a glowing combustion state to a flaming combustion 
state in scaling with do: 


In each of Figure 7 to Figure 14 the curves of the four 
irradiance levels effectively collapse to a single curve. The data 
points coalesce especially well in the beginning and the ending 
Part of the burn history in each’ curve, thus supporting the 
reasonable assumptions made in formulating the time stretching 
parameter, G,- In summary, the utilization of a fitting parameter, 
qo? and of the time stretching parameter, Gi» was effective in 
collapsing the data at different irradiances to a single curve. 
Since these parameters were also based on a physical principle, 
although not in a rigorous fashion, we can interpolate and extrapo- 
late to a certain extent from the scaled burnrate database as a 
molvou1 on Of .the net surface yneat flux .and.burn .nistory... More 
specifically, by rearranging the scaling relationships in Equations 
(3) to (9) and noting the results in Figures 7 to Figure 14, it is 
still possible to define a heat of pyrolysis as a function of 
scaled time and define a heat of combustion also as a function of 


scaled time as shown in Equations (10) and (11) below. 


x 
Peaumou pyrolysis -is function of ,T, only: 


q ° 
aan eo 2 On- e eerg  ee (10) 
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for polyurethane: 


q 20 kW/m? 


e) 


4 
i) 


LOUPK 


This approach gives us more confidence to interpolate and extrapo- 
late from the cone calorimeter burning data at different 
irradiances to a time changing burn rate of a furniture surface 
element. More burnrate databases are being incorporated into the 
computer files in order to more severely test the scaling relation- 
Slerpisi 


oss FLAME STRUCTURE MODELING IMPROVEMENTS 


Both the flame shape and the flame properties were changed in 
the flame structure module as a result of new information. ‘The 
changes in this module did not affect the other main modules be- 
cause of the geometries chosen for the furniture fire model. Each 
furniture cushion panel can have several burning surface elements, 
each having its own fuel loss and heat release rates as a function 
of burn history, and the panel can also have a few burning polygons 
which represent the flame spreading process. The cushion panel is 
a quadrilateral plane which can have any position or orientation in 
a three dimensional space. Four such cushion panels are combined 
to construct a generic chair. “A rigorous three-dimenslonal -veoeor 
analysis was utilized to represent the geometries of the furniture, 
flame, and the room for the reasons of flexibility, compactness, 
and-efficiency. Thus it was quite easy to “replace ‘the *conreas 
frustrum flame shape used in the demonstrator version by a polygon- 
based parallelpiped shape shown in Figure 15. The new shape was 
decided after our recent viewing of the VCR recording of the Pura 
CUrewMOCKUD [Timer tests atyNBs. 


Each flame shape is sub-divided into surface elements. Each 
Surface element on the chair, the flame, or the room is charac- 
terized by its three-dimensional position vector, the element 


normal vector, and the element area. These characterizations are 
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Vertical Burning Polygon —>>» 


Horizontal Burning me) 


Figure 15. Flame Shape is a Polygonal-Based 
Parallelepiped. 
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used in the thermal radiation analysis to calculate the viewlagveee 
and the Hottel exchange areas between all surface elements. The 
mean beam length approximations utilized in the Hottel exchange 
areas required computations of the flame volume, the flame surface 
area, and the burning polygon centroids, areas, and circumferences. 
The three dimensional vector analysis is used to obtain these 
geometric parameters. The mean beam length approximations also 
implied the use of the flame properties of soot and gaseous vex. aie. 
tion coefficients and of the flame temperature, which are assumed 
spatially uniform inside a flame shape. This is consistent with 
the formulae presented for the flame properties by Bard and Pagni 
(Reference 4). The flame height formulae are presented in 
Equations (12) tof (iS):<pelow. 


Cetegen, et. al., formula for averaged visible pool fire height is, 


(9*) seg 
ee Qp ; Qn 4 
Liss ALD (12) 
p 
ele % 
where 
pe Vh285/2 
Sie (Q,/Le,C,T.8 ECO 4) 
D = hydraulic diameter 
wa = total heat release rate. 


Hasemi's formula for averaged visible wall fire height Ls? 


4Y/5 
x ¥ 
6.0 (Q,) sy ean 
train 3 (13) 
: eee cS % 
0 by : Qy >a 
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where 


i If 2s Ave 
Q. (Q,/[o.C Tog Resbenhoslp 
ZN = pyrolysis zone length 
oF = pyrolysis zone width 


Interpolation for merged flame: 


ee 9 ae? (14) 
fi fo Qa fr Q, 
where 
Q, = heat release rate from horizontal surfaces 
p 
Qa = heat release rate from inclined surfaces 
W 


with a solid flame height: 


Z = 0.4 Ze CLS) 


The Cetegen, Zuboski, and Kubota's formula (Reference 8) for the 
averaged visible pool fire height was used for the conical frustrum 
flame shape in the Demonstrator version. This formula is also used 
in the upgraded Flame Structure Module along with Hasemi's recent 
formula (Reference 9) for the averaged visible wall fire height. 

It is noted the pool fire height is sensitive to the total heat 
release rate and to the hydraulic diameter of the flame polygonal 
horizontal base, while the wall fire height is sensitive to the 
total heat release rate and to the length and width of the inclined 
burning polygon. The effect of the flame radiative fraction on the 
flame height is considered minor’ and is. not? ineluded at this stage. 
However, if both the seat and the back rest on the furniture are 
burning, an interpolation between the pool and the wall fire 


Zo 


heights is required and is given by Equation (14). The linear 
interpolation weights the pool fire height with the heat release 
rate from horizontal surfaces and weights the wall fire height with 
the heat release rate from inclined surfaces. The resulting 
visible height of the furniture flame shape is further reduced by 
Sixty percent to obtain the 'solid' flame height defined in the 


above two references. 


2.4 THERMAL RADIATION ANALYZER MODULE 


Improvements in the flame properties module were mainly in 
the calculation of the soot and gas extinction coefficient. angeine 
corresponding flame emissivity. Using the results from Bard and 
Pagni (Reference 4) for radiation from soot and gas for any size 
pool fire, the inputs needed are the maximum soot volume fraction, 
the partial pressures of water and carbon dioxide gases, the mean 
beam length, and the flame temperature to calculate the extinction 
coefficient as given by Equations (16) to (23) below. 


Bard and Pagni's scaling relationship for flame properies of pool 


fires: 
qd = €,0 rt (16) 
ie f 
Daca: age CT#2 
—€ = 1 - exp(-kL) (18) 
K = 367 boa.) f/h (19) 
AT p = 0.40 cm°K (20) 
oa 4 ia’ (21) 
U3 y a es ee, ; Kise Se V3. 
¢; ite = ie KL 2 Q .43 
max 
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For polyurethane, they obtained: 


Ft 
i] 


“4 0.89 ppm 22") 


J 
i] 


f teu (230) 


The soot extinction coefficient is a function of a particular mean 
wavelength and of the soot volume fraction. The mean wavelength is 
inversely proportional to the gas radiative temperature. The soot 
volume fraction reaches a saturation level for a large enough pool 
fire. For these reasons the model for the soot and gas extinction 
coefficient for the upper gas layer in the FAST model was also 
examined. The flame temperature of the pool/wall fire was assumed 
to be at 1200° Kelvin and the maximum soot volume fraction was 
taken to be 0.89 ppm. as corresponding to polyurethane foam 
materials. | 


Once again a main advantage of the furniture fire model is 
the flexibility in being able to completely revise any particular 
main module listed in Appendix A without affecting the other 
remaining modules. The freedom to choose any particular flame 
Shape and to use a realistic flame property data is a result of 
relying on efficient computations of viewfactors between all sur- 
face elements as a basis for evaluating the Hottel's exchange areas 
between the furniture, room, flame, and the gas layers. As the 
Surface elements are made smaller and as the flame approaches the 
optically thick limit, the radiative heat transfer approaches the 
theoretical solution on a sufficiently large computing facilities. 
Yet the furniture fire model is capable of running on a advanced 
personal computer without the surface elements becoming too big for 
a significant loss of accuracy. To allow this flexibility in 
adaptation to a computer facility the user can provide inputs to 
the geometrics of the model to control, (a) surface element size, 
(b) the resolution and the number of burning polygons, and in a 
later version (c) constructions of furniture shapes and multiple 
flame shapes. 
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SECTION 3 
HEMFAST CODING DEVELOPMENT AND RESULTS 


The emphasis during this contract period was on the inclusion 
of the furniture fire model into the model for the Transport of 
Fire, Smoke, and Toxic Gases (FAST) program. The flowchart for the 
main modules of HEMFAST are located in Appendix A. All modules 
shown were developed at UDRI except for the FAST model. Extensive 
modification of the thermal radiation and the fire source modelings 
in FAST were performed to accomplish this task. The flowchart of 
each module of the furniture fire model is also included in 
Appendix A. Each furniture fire main module has as input the 
processed bench scale data. . These data are described in tnhestor. 
lowing paragraphs. The merged program is now identified as the 


Heat and Mass Transport of Furniture Fire, Smoke, and Toxic Gases. 


A magnetic tape copy of the HEMFAST program and its input 
file was submitted to NBS. A hard copy of the sample run was also 
provided. The computer programmer's comments and suggestions on 
executing the HEMFAST program on a CDC Cyber, such as the one at 
NBS, are located in Appendix B. 


See DISCUSSION OF FAST MODEL 


The FAST model is a room model developed and presently used 
at NBS Center of Fire Research. The model employs mass flow zone 
modeling of smoke, soot, and molecular constituent movements in a 
multi-room configuration. The model presently has a diathermious, 
two-layer, and gray gas formulation of the radiative heat flux 
taken from Engel and Howell (Reference 10). The 3-D Hottel zoning 
formulation presently used in the furniture fire model was easily 
extended to the room environment by appropriate geometrical 
modifications in the computer program. Thus it was ossible to 
calculate the radiative heat flux inputs for the four room bound- 
aries: ceiling, upper wall, lower wall, and the floor, to the time 
dependent heat conduction solution of the surface temperature. 


at 


The FAST model presently uses as input, data from a gas 
burner or other fire source measured in the laboratory. Therefore, 
the model validation of FAST is limited to a few particular 
laboratory experiments until a self-consistent fire source model is 
developed. The incorporation of the furniture model into FAST was 
a way of providing such a self-consistent fire source. 


Version 17 of FAST was sent to UDRI in December 1985. This 
version was modified for incorporating the furniture fire model. 
Examination of the FAST code revealed that only seven subroutines 
of the FAST model required modification. Six of the subroutines 
are called by the seventh subroutine, named SOURCE. There is a 
block of approximately 30 lines of code in SOURCE that pertain to 
the radiative heat flux and the furniture fire models alluded to 
earlier. Instead of replacing this block of lines in the SOURCE 
Subroutine, we developed an optional bypass routine with a similar 
number of lines of code that will call the furniture fire model and 
call the extended 3-D radiative heat flux model to take into ac- 
count the upper gas layer, the fire source, and the geometrical 
arrangement of the furniture and room in one set of calculations. 
The furniture fire model provides the total burnrate and the mass 
fractions of soot and of other species constituents as an input to 
FAST for computations of the upper gas layer. 


Possible interaction problems which included consistencies in 
units, coordinate systems, and the time-stepping scheme were 
resolved. Time-stepping was controlled in such a way to avoid 
instabilities, since the present furniture fire model relies on an 
explicit time integration scheme. 


Extensive, but straightforward, reevaluation of Hottel's 
exchange areas were needed to incorporate the upper and lower gas 
volumes and to incorporate the room boundaries. Since all surface 
elements in the room were included, this gave rise to radiation 
blocking between the chair and the room. With the extra gas layers 
and the extra surface elements defined, the smoothing of Hottel's 
exchange areas to preserve energy conservation was severely tested. 


Its algorithm also was revised. 
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Bee HEMFAST DATA INPUT DESCRIPTION 


There are two input data files used by the HEMFAST model. 
The input data file is essentially in the format’ of Vers ior). 
the Transport of Fire, Smoke, and Toxic Gases (FAST) program and a 
description of the data file is given in the User's Guide (11). 
For the sample run, a rectangular room with a single window-opening 
to the outside was ‘geometrically simulated. = The \conductivesane 
emissivity properties in SI units were provided for the ceiling, 
wall, and floor. When the computer: index variable,\:(LFBT is ;Segeva 
5, the type of fire simulated is the furniture fire. For LFBT 
Other than five, the HEMFAST program reduces to the FAST program. 
Thus, the modularity of fire types is retained, although the 
program code may need further restructuring and inner documentation 
to make it more convenient to the user. Since the furniture fire 
model is a self-controllable fire source the value of LFMAX in 
input data file need not be greater than one. If the value of 
LFMAX is one, it can be used to define the time interval of fire 
growth in seconds and to initialize the mass fractions of the 
product gases of fuel mass combustion. These initial values of 
mass fractions were held constant throughout the fire growth if the 
corresponding bench scale data of the molecular constituents were 
not available as a function of time. The initial values of mass 
oD CO,, Coe HO, and soot were taken to be that if a 
generic polyurethane (PU) foam. 


fractions of N 


The second input data file contains bench scale fire 
properties of a furniture material and the specifications for a 
four cushion chair or couch geometry. The furniture material 
selected for the test case was a 50mm thick, light 
olefin/California PU foam material (3). The bench scale data 
consisted of: 


a) -. burn rate vs time for o5kW/m* irradiance of the 
b) - mass loss rate vs time cone calorimeter (Reference 3) 
c) - ignition temperature (711°K) (Reference 12) 


d) - thermal thickness (0.72 KJ/°Km*) (Reference 12) 
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e) - horizontal flame spreading constant of the material 
(Reference 12) 


f) - maximum soot volume fraction (0.89 PPM) (Reference 4) 
g) - flame temperature (1200°K) (Reference 4) 


The furniture geometry simulated was the four-cushion furniture 


mockup similar to that shown in Krasny and Babrauskas (3). A 
Simulated cushion was considered to be a 400 mm length panel with 
four surface elements. Four of the cushion panels were connected 


and arranged into a chair for the demonstration computer run. The 
user can control the chair geometry with the input data file. The 
last entries in the data file refer to additional input in the 
Spark ignition at time "zero" and its location on the righ:-rear 
surface element of the seat panel. 


3/3 HEMFAST DEMONSTRATION RESULTS 


As indicated in the FAST User Guide, the first output data is 
a summary of input data. When the fire type is the furniture fire, 
the FAST Output was changed so that the "Time Step" is now control- 
led by the user and "Number of Fire Intervals", "Effective Heat of 
Combustion", "FMASS", "FHIGH", and "FTIME" output parameters are 
ignored. In later editions of the code these types of outputs will 
be suppressed. At "Time = 0.0 seconds" the output variables of the 
SoLiculations are initialized: Note that concentrations of N, and 
HO are initialized. The amount of water vapor assumed to exist in 
the lower gas layer corresponds to 50% relative humidity at a room 
temperature of. 300°K. The nexi.setr of eoutput is relatedyto. the 
description of the furniture surface elements located relative to 
the room origin (the right-rear corner of the room as shown in 
Figure 16). Each surface element's mid-point vector, area, and the 
normal vector are needed as inputs to the thermal radiation model 


of the furniture fire model. 


At the first time increment, (10 seconds) the furniture fire 
sub-model computes a pool flame (centered at 2.2, 2.05, 0.5 meters) 
on the first surface element of the furniture with a base radius of 


a5 
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Fig. 16. Arrangment of Surface Elements 
on Chair and Room. 
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0.01 m, flame radiative height of 0.001 m, and a burnrate of 

0.001 kW (see Figure 17). The result of the radiative heat trans- 
fer analysis provides interesting and reasonable features. 
Included in the printout are the three isothermal gas volumes, the 
pool flame, the upper gas layer, and the lower gas layer. Also 
included are all the surface elements of the furniture and of the 
room. The first volume element, the pool fire volume, has the 
following result: zero scattering, a small total fire area, small 
gray-body emissivity of 0.0027, a black-body power of 11.76 kW/m, 
Slimereodiance of 0.472 kW/m* from the upper and lower gas volumes 
and all the surface elements, a radiance of 117.6 kW/m", and a net 
radiative power loss of 0.00021 kW which was 21% of the burnrate. 
The small upper gas layer has a net radiative power gain of 0.11 x 
10. ‘kW. The lower gas layer on the other hand has a large gray- 
body emissivity of 0.143 due to the large mean beam length and the 
presence of 50% relative humidity. It also has a large net radia- 
tive power gain (0.00056 kW) compared to the other elements. The 
only element with a higher net power gain (0.00058 kW) is Element 
ieeetne ceiling. 


The first surface element on the furniture contained the 
relatively small fire base, which resulted in a slightly higher 
Overall surface emissivity of 0.8014. However, its overall 
radiance of 0.541 kw/m° was higher than the irradiance of 0.461 
kw/m-, resulting in a net radiative power loss of 0.032 kW. This 
LSamuchy greater than »themet power) loss.of theepoolyfire. This is 
not too surprising considering the small gray-body emissivity of 
the pancake shaped pool fire. To sustain the ignited surface 
element base at 711°K and a radiance of 0.0032 kW requires glowing 
combustion of the fabric/foam material. All the other surface 
elements should be absorbing net radiative power in relation to 
their view factors from the pool fire. This is generally the case, 
except for three insignificant surface elements. These exceptions 
are the result of either roundoff errors or an error tolerance of 
the Gauss-Siedel matrix solution not being small enough. 
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Fig. disso HEMPASTY Prediction of puri eee an 
versus Time for the Sample Chair Fire. 
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These overall results of the radiative heat transfer are 
unexpected of large pool fires, but they are essential to the early 
development of the furniture fire. For the rest of the simulation 
at time = 10 seconds, not much was changing in the other variables. 
All the temperatures were still at approximately 300°K (see Figure 
18), except for the first furniture surface element midpoint, which 
was at ignition temperature of 711°K. The upper layer species 
concentrations of COo, COT HO, OD and CT started at very low 
values (See Figure 19). 


The HEMFAST results began to change in character by time = 40 
feconds. =< The pool fire was then burning at a rate of 2.1 KW, as 
Shown in Figure 17, the fire base radius was 0.09 m, and the flame 
magiavive height was 0.05 m. The net radiative fraction, of the 
pool fire burnrate has increased to 25% and the flame emissivity 
has increased to 0.05. The flame has a net radiative power loss of 
0.51 KW, which at this time is greater than the net radiative power 
loss of 0.3 kW of the ignited surface element. The upper and lower 
gas layers were experiencing net radiative power loss instead of 
net power gain. This was due largely to the temperature of the 
upper and lower gas layers being 10°C greater than the temperature 
of the ceiling, walls, and floor (see Figure 18). The upper gas 
layer was 0.5 m deep with a temperature of 310.9°K. The parts per 
million concentration of CO., CO, and H50 werew~isc, 6.0, and 780, 
respectively (see Figure 14). The soot opacity was 6 x 103 per 
meter. The armrest surface element adjacent to the burning seat 
element has its temperature at 324°K and the adjacent backrest 
surface element was at 321°K. 


At time = 50 seconds, with a burnrate of 3.93 kW, the pool 
fire has spread over the adjacent surface elements on the armrest 
and backrest. As a result, the pool fire radiative height began to 
rapidly increase to 0.15 meters and the radiative fraction in- 
ereased to 45% over the next 20 seconds. At 50 seconds, the pool 
fire has become the dominate radiation emitter at 1.23 kW and most 
burning surface elements still were experiencing net radiative 
power loss. At this point it is uncertain whether the burning 
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surface elements were sustained by glowing combustion (because 
convective heating was not yet directly calculated), nor was it 
certain the burning surface elements will remain at the ignition 
temperature of 711°K. 


Indeed, at 70 seconds a:small problem occurred when the 
temperature prediction of some virgin surface elements dropped 
Significantly below 300°K. This may be related to: (1) a 
stability problem, thus requiring smaller time step; (2) not con- 
sidering a wall fire as separate from the pool fire; or (3) more 
explicit use of convective heat transfer calculations may be 
needed. A numerical stabilizing influence occurs when the surface 
temperature was predicted to rise above the ignition temperature. 
This indicates spontaneous thermal ignition, which would happen 
with a more stable surface temperature prediction. The surface was 
then set at the ignition temperature, 711°K, and the burning rate 
initialized for that element. This, in fact, happened at time = 80 
seconds on the left armrest and then at. time = 90 seconds on the 
remaining surface elements on the backrest. 


At time = 80 seconds, the burnrate was 20 kW, the whole seat 
and parts of the side and backrest were burning, the flame radia- 
tion height was 0.15 m, and the radiation fraction became 51% of 
the burnrate. With the flame emissivity approaching 0.18 with a 
radiative net power loss of 10.3 kW, some burning surface elements 
began to switch to net radiative power gain. This in turn enhanced 
the burning rate of the surface elements. 


At 90 seconds, the burnrate was 52 kW, the flame radiation 
height became 0.3 m, and the whole mockup was on fire. The flame 
emissivity continued to increase up to 0.28 with radiative net 
power loss of 23.6 kW. At this point in time the fire was ex= 
periencing its greatest growth rate, that. is, all the furnimae 
Surface elements had ignited and were burning. The flame spreading 
process has nearly finished and the continued fire growth is 
resulting from the increasing flame radiative heat feedback to the 
burning surfaces. Eventually, the fire growth reaches a limit 
dictated either by the flame soot volume fraction reaching its 
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maximum value for the given flame base radius or by the burnout of 
Some surface elements. 


Other calculated features of the furniture fire growth at 90 
seconds are the following. The upper gas layer has reached a 
limiting depth of 0.7 m considering the open window. The tempera- 
ture of the upper gas layer has reached 346.5°K, while the lower 
gas layer has reached 331.3°K mainly because of the radiant absorp- 
tion by the water vapor. The ceiling, upper wall, lower wall, and 
floor temperatures were calculated as 309.4, 306.2, 306.0 and 
300.5°K, respectively. The upper layer species concentration of 
CO, CO, and H,0 were 29707" 129) and416960U ppm, sréspect ively: ithe 
Soot opacity was 0.227 per meter. 


At 110 seconds, the simulated furniture burnrate finally 
reaches a peak of 109 kW. This gives a flame radiance height of 
0.5 m and a flame gray emissivity of 0.326. The net radiant power 
loss was 40 kW for the flame, 2.6 kW for the upper has layer, and 
-1.1 kW for the lower gas layer. It is interesting that the upper 
gas layer consistently had a net radiant energy loss to the lower 
gas layers and to all the walls despite the radiant heat gain from 
the furniture fire. Yet the upper gas layer was still heated by 
the convective heat portion of the furniture burning rate. This 
process gave the result of 387°K for the temperature of the upper 
weomrayer With ts depth stillv*at O77 m-*etne calculated ‘tempera- 
tures for the lower gas layer, ceiling, upper wall, lower wall, and 
Pieseroor were’ 33829°K, 315. 21°K, “310 -6°K, °309,0°K/ and 307. 1°R, 
respectively. The upper layer species concentration for cO., co, 
and H.0 reached their peak values of 4997.0, 217.0, and 19640 ppm, 
respectively. There were a couple of furniture surface elements 
that burnt out, thus indicating a rapid downward trend of the 
furniture burning rate. By the time 130 seconds, the mockup was 


DUrNe oat: 
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SECTION 4 
CONCLUSIONS 


HEMFAST is a computer program which simulates the furniture 
fire and the toxic gas growth in a room and its spread to the other 
rooms. This was achieved by incorporating the furniture fire model 
into the FAST model. The furniture fire model was designed to be 
self-controllable; in,»simulating: a:  fire,source utilizing. aseinguTs 
only the bench scale data of material burning rate, flame spread- 
ing, thermal ignition, and the flame properties of the following: 
soot volume fraction, COs, H,0 mass fractions of fuel weight, and 
gas temperatures. The Froude modeling of flame shape and of 
uniform flame properties is utilized following McCaffrey (13) and 
Orloff.and de Ris (6)...To. obtain the self-controllabil ityeeamee 
furniture fire model, and also in HEMFAST, a practical 3-D thermal 
radiation model was used (14) in calculating the heat feedback from 
the flames to the flame front in the determination of the flame 
Spread rates, and to the flame base in the determination of the 
local material burning rates. In the HEMFAST program, tne Unerman 
radiation model (14) was extended to include the upper and lower 
gas layers and the room enclosure. Since the Hottel's zoning 
method was used in the thermal radiation model it was possible to 
consider, surface emissivities less than unity on all surfaces weang 
to consider three isothermal gray gas volumes which are the high- 
temperature flaming plume, the sooty upper gas layer, and the lower 
gas layer at 50% relative humidity. 


Some salient features were noted in the simulation of a 
sample fire growth on a furniture mockup. Immediately after the 
seat was ignited, the pool fire was inferred to be sustained by 
glowing combustion within the burning material while the flame 
front was spreading at its slowest rate. This was shown by the 
relatively high net radiative power loss calculated for the burning 
material. The simulated furniture fire tended to grow at a low 
exponential rate for a slow but increasing flame spread rate from 
the octogonal flame base. 
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As the flame base reached the right armrest and the right 
Side of the backrest, a new phase in the fire growth rate 
developed. This new phase is caused by the additional burnrate of 
the side walls and its radiative feedback to the burning seat. The 
increased burning rate led to a higher flame radiative height and 
also to a higher soot volume fraction. The final phase in the fire 
growth rate involved the spontaneous thermal ignition of the 
remaining surface elements. This phenomena has also been observed 
in VCR observations of mockup burning. During the first two phases 
fener: ire vrowthn, i.e. ,° the seat and side wall flame spreading 
many of the burning surface elements had a net radiative power 
loss. Only after the spontaneous thermal ignition of several 
Surface elements was the radiative feedback sufficient to the 
result in net radiative gain for these burning surface elements. 
With all the elements burning at this time, the fire growth was at 
its greatest. However, the peak in the furniture fire burnrate in 
this case was limited by the surface elements burnout as opposed to 
being limited by the maximum soot volume fraction for a given flame 
base size. 


The upper and lower gas layers eventually rose 87 and 39°C 
respectively above the initial room temperature while the upper 
5» COs H50 eventually reached 4997, 217, 
and 19640 ppm respectively at peak burning rate. Throughout most 


fayer constituents for CO 


of the fire growth, the upper layer had a net radiative heat loss 
to the lower gas layer and wall despite the heat gain from the 
furniture fire. Whether the upper gas layer will have a net radia- 
tive heat gain or loss depends on its geometry, soot volume 
fraction, and gas temperature in relation to other radiative 


sources. 
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SECTION 5 
RECOMMENDATIONS 


Although the HEMFAST program currently predicts the qualita- . 
tive behavior of the furniture fire growth in a room, it will need 
some specific revisions for quantitative comparisons with the data, 
especially the data by Krasny and Babrauska (3). Also, the current 
version of HEMFAST is not user friendly; is not thoroughly 
debugged; does not yet incorporate all the features; and is not 
thoroughly optimized. It would be useful to somehow derive the 
flame properties of soot volume fraction, CO., and HO mass frac- 
tions, and the gas temperatures within the flame from the bench 
scale data. Currently, these flame properties are from the data 
obtained by Pagni and Bard (4) for a PU foam. The bench scale data 
available on the horizontal and lateral flame spreading on furni- 
ture materials is very limited. At the present some data on the 
flame rate can only be inferred from the VCR observations of furni- 
ture mockup burnings (3). Also, for some cases, the data for the 
temperature of the burning surface is nonexistent. The HEMFAST 
model can be improved in the following areas: (1) the prediction 
of the furniture surface temperatures by using more stable numeri- 
cal technique, such as the one used in the FAST model for the room 
wall temperatures; (2) insert an additional flame source by includ- 
ing a wall fire algorithm of existing integral models of wall fire; 
(3) perform a more in-depth processing of the bench-scale data on 
the burning rates of furniture materials as measured by the 
Calorimeter to infer the burning surface temperature and scaling 
factor due to the different irradiances from the cone calorimeter; 
(4) do a thorough re-examination of the Hottel's exchange area 
approximations in the thermal radiation model especially regarding 
the conclusions about the importance of glowing combustion in the 
early phases of the furniture fire growth, and (5) the computer 
programming of HEMFAST can be improved to make it more user 
friendly, more modular, and more extensively documented. An area 
of documentation needed is in geometrical representation of the 


furniture and flame shapes by surface elements and of the flame 
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Spreading regions by surface polygons. The role of three- 
dimensional vector analysis in deriving geometrical properties such 
as flame volume and areas as well as flame base size and area 
should be included in the documentation. 
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Input: Local Ip, I¢ 


Surface Temperature Model 
for Virgin Material, TS and e. 


Surface Element is 
Ignited and cs is 


given 


Pyrolysis Boundary on 
Surface Element? 


tocal. TT) > T, ? 
rae He | 


Output: Local e. and T> 


Flame Spread Velocity Model 


Vv ate 
p Pp 


as a Function of Time 


Pyrolysis Boundary Model 


r= f v_dt 
P P 


Output: local e., TL, and £ 


as a Function of Time 


Local Surface Temperature Ignition, and Flame Spread Predictors 


A-3 


Input: Local da and os 


Model of T , ee, andm 
Ss s Dp 


for Burning Material 


Burning Material Thickness 
Model, 6, = [mdt/o, 


| Set qe 70 and e,70 | 


Burnout Models: 
Surface Elements 
Removed and Burnout 
Burning Rate Formula Boundary Redefined 
e > 


for qo Eb 


Output: Local e,, Tg, qe and Ss 
as a Function of 1e 


Local Burn Rate, Fuel Consumption, and Convective Heat Flux Predictors 


Burning Area Definitions 
and Surface Orientations 


Gp >a 1 
Isolated Tall Tsolated 
Pool Fires Wall Pires 
BN Scrit 


Isolated Short Tsolated 
Pool Fires Corner Fires 


Plame Merging 
Model 


Overall 
Plame Structures 


Flame Structure Model 


Flame and Body 
Structure 


3-0 Confisuration 
Factor with Blocking 


Hottel's Exchange 
Areas Computation 


Room Exchange Areas 
via Energy Conservation 


Radiation Equations Solution 
via Gauss-Seidel Iteration 


ov and Qi Output 


Thermal Radiation Heat Transfer Analyzer Module 
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APENDIX B 


PROGRAMMER'S COMMENTS 


PROGRAMMER'S COMMENTS: 


The tape is written at 1600 BPI, 

Block Size is 5120 

Record size is 80 characters, fixed length 
No carriage control 


The tape contains three files: 


1. The HEMFAST Program 
2. Data file A (FAST.DAT) 
3. Data file B (FURNITURE. DAT) 


Data file A contains the data required by Version 17 of the 
NBS FAST model 


Data file B contains data required by the UD Furniture Fire 
model 


Data file A is read from FORTRAN Unit 001 
Data file B is read from FORTRAN Unit 005 


The output is written to FORTRAN Unit 006 (a sample output 
is enclosed) 


you will probably need to add two or three OPEN statements 
at the beginning of the program to take care of this. 


APPENDIX C 


PRINTOUT OF HEMFAST RESULTS FOR 
TEST CASE FURNITURE FIRE AT 
SIMULATION TIME, 60 SECONDS 


BUILDING XYZ 


TOTAL COMPARTMENTS = 
MAXIMUM OPENINGS PER PAIR = 


FLOOR PLAN 
WIDTH 4.6 
DEPTH 4.3 
HEIGHT 2.4 
AREA 19.8 
VOLUME 47.35 
CEILING 2.4 
FLOOR 0.0 
CONNECTIONS 
1 ( 1) BW= 0.0 pi 
HH= 0.0 2. 
HL= 0.0 QO. 
HHP = 0.0 2. 
HLP= 0.0 0. 
CEILING 
COND = 1. BOOE-04 
SPHT = 9. OOOE-01 
DNSTY= 7. 9OOE+02 
THICK= 1. 600E-02 
EMISS= 9. OOOE-01 
FLOOR 
COND = 1. 7OOE-03 
SPHT = 1. OOOE+00 
DNSTY= 2. 200E+03 
THICK= 1. SOOE-O1 
EMISS= 9. COOE-O1 
UPPER WALL 
COND = 1. BOOE-04 
SPHT = 9. OOOE-01 
DNSTY= 7. 9OOE+02 
THICK= 1. 600E-02 
EMISS= 9. OOOE-01 
LOWER WALL 
COND = 1. BOOE-04 
SPHT = 9. OOOE-O1 
DNSTY= 7. 9700E+02 
THICK= 1. 600E-02 
EMISS= 9. COOE-O1 


FIRE ROOM NUMBER IS 1 

TIME STEP IS 10.00 SECONDS 
PRINT EVERY 1 TIME STEPS 
NUMBER OF FIRE INTERVALS = 


oOreoreW 


TOTAL TIME INTERVAL = 18 
FIRE SOURCE = 1 
FIRE TYPE = FURNTRE FIRE 


INITIAL FUEL TEMPERATURE (K) = 300. 
AMBIENT AIR TEMPERATURE (K) = 300. 
AMBIENT REFERENCE PRESSURE (kPa) = 101.30 
EFFECTIVE HEAT OF COMBUSTION (KUJU/KG) = 18000. 
FMASS= ©. OO0E+00 0. 0OE+00 
FHIGH= 0. 0O0E+00 0. OOE+00 
N2= fal 0. OOE+00 
CO2= 2.2 0. OOE+00 
CO= 6.00E-O2 OQ. O0OE+00 
H20= 0. 67 0. OGE+00 
OD= 2. 00E-02 0. O0E+06 
CT= 1.0 0. QOE+0C6 
FTIME= 10. 
-ONFG FILE = 
FURNITURE BURNRATE IS 12. 36035 KW AT TIME 60. 00000 SEC 


SMFAC < O.S IN RMTERM SUBROUTINE MEANS A 0. 3953227 
tELATIVE CHANGE IN THE SUM OF SURFACE-GAS EXCHANGE AREAS 
PROVIDED THE SUM OF GAS-GAS EXCHANGE AREAS ARE CORRECT; 
THE SUMS WERE FOR THE GAS VOLUME #¢ 1 


RESULT OF RADIATIVE HEAT TRANSFER ANALYSIS FOR THE NEXT TIME INCREMENT AT FIRE 
SOURCE ROOM 

THE VAPOR PRESSURE FOR CO2 FOR LOWER AND UPPER GASES AND FOR H20 FOR LOWER AND 
UPPER GASES ARE 0. 0000000E+00 2.3781539E-04 1.5071809E-02 1. 0837992E-02 


ITERATION FLAME. AREA RAD. FRACTION FLAME. HEIGHT FLAME. RADIUS 
4 0. 2683476 -0. 3345792 0. 1523599 0. 1289839 
ELE SCAT/COSAV ELEM. AREA EMISSIVITY BLCKBDYPOW IRRADIANCE RADIANCE NETPOWER (KW: 
1 0. OOOOE+00 O. 2683E+00 0. 1541E+00 0. 1176E+03 0. 2640E+01 0. 1176E+03-0. 4136E+0) 
2 0. OOOOE+00 0. 4924E+02 0. 8132E-01 0. 5627E+00 0. 5357E+00 0. 5627E+00-0. 1081E+0¢ 
3 0. OOOOE+00 0. 7416E+02 O. 1271E+00 0. 5648E+00 0. SS5S0E+00 0. 5S648E+00-0. 2811E+0¢ 
1 0. 1O000E+01 0. 4000E-01 0. 1Q00E+01 0. 1451E+02 0.1170E+02 0. 1451E+02-0. 1125E+0¢ 
2 0. 7223E+00 0. 4000E-O01 O. 8O00E+00 O. 1860E+01 O. B381E+01 0. 2164E+0!1 O. 4865E-0: 
3 0. 7239E+00 0. 4000E-01 0. 8000E+00 0. 1954E+01 0. SB97E+01 O. 2342E+01 O. 6218E-0: 
4 0. 4525E+00 0. 4000E-01 O. SO00E+00 O. 1348E+01 0. 371I9E+01 O. 1822E+01 O. 7589E-0: 
5 0. 5550E+00 0. 4000E-01 O. 1OO0E+O!1 O. 1451E+02 0. 6337E+01 O. 1451E+02-0. 3271E+0¢ 
6&6 0. 4540E+00 0. 4000E-01 0. 8000E+00 0. 1777E+01 0. 4268E+01 0. 2275E+01 0. 7973E-0: 
7 0. 2262E-01 0. 4000E-01 O. 1000E+01 O. 1451E+02 0. 7O15SE+01 O. 1451E+02-0. 3000E+0¢ 
8 0. 6487E+00 0. 4000E-01 O. B8000E+00 0. 2218E+01 0. SVO4E+O01 O. 2955E+01 0. 1179E+0¢ 
9 0. 1O00E+01 0. 4000E-01 O. 1O00E+01 O. 1451E+02 0. 1214E+02 0. 1451E+02-0. 9495E-0: 
. 6942E+00 0. 4000E-01 0. BO00E+00 O. 1691E+01 0. S586E+01 0. 2470E+01 O. 1246E+0( 
. S496E+00 0. 4000E-01 O. 1O00E+01 O. 1451E+02 0. S406E+01 0. 1451E+02-0. 3643E+0¢ 
. 4553E+00 0. 4000E-01 0. SOOOE+00 O. 1362E+01 0. SB90E+01 0. 1867E+01 0. BOB9E-0: 
. F262E+00 0. 4000E-01 O. 8OO0E+00 0. 7332E+01 0. 6400E+01 0. 7145E+01-0. 2983E-0: 
. BOFVOE+00 0. 4000E-01 0. BOO0E+00 O. SBBSOE+01 O. S631E+01 O. SB30E+01-0. 7952E-0: 
. 7512E+O00 0. 4000E-01 0. BO00E+00 0. 6992E+01 0. S532E+01 0. &6700E+01-0. 4669E-0: 
. 6751E+00 0. 4000E-01 0. SOOOE+00 0. 6262E+01 0. SO9SE+01 0. S5030E+01-0. 3722E-0: 
. 99O03E+00 0. 1978E+02 0. SOOOE+00 0. 4718E+00 0. 5323E+00 0. 4839E+00 0. 9580E+0( 
. 7B18E+00 0. 2338E+01 0. BOOOE+00 0. 4688E+00 O. 5S253E+00 0. 4801E+00 O. 1057E+0( 
. 7968E+00 0. 2501E+01 0. BOOOE+00 0. 4688E+00 0. 5471E+00 0. 4845E+00 0. 1567E+0¢ 
. OOOOE+00 0. 2338E+01 0. BOODE+00 0. 4688E+00 O. 5S326E+00 0. 4816E+00 0. 1194E+0( 
. BOOSE+00 O. 2501E+01 0. BOOD0E+00 0. 4688E+00 0. S029E+00 0. 4756E+00 0. 6814E-0: 
. 9829E+00 0. 79B82E+01 0. BOOOE+00 0. 4683E+00 0. 5143E+00 0. 4775E+00 0. 2936E+0( 
. 9796E+00 0. B8S39E+01 0. BOOOE+00 0. 4683E+00 O. 5682E+00 0. 4883E+00 0. 6827E+0! 
. OOOOE+00 O. 79B2E+01 0. BOOOE+00 0. 4683E+00 0. 4972E+00 0. 4741E+00 0. 1844E+C0 
.9796E+00 0. BS39E+01 0. BOOOE+00 0. 4683E+00 0. 5157E+00 0. 4778E+00 0 3240E+C 
. 8989E+00 0. 1978E+02 0 B000E+00 0. 4597E+00 0. 6074E+00 0. 4893E+00 0. 2336E+C 
IN KW ABSORBED BY ALL SURFACES I 4. 4981 


TIME = 60.0 SECONDS. 
U. TEMP. 322.8 
CATEMP.. 320. 0 
U. VOLUM 11.9 
U. DEPTH 0.6 
CEs TEMP 303. 1 
UW. TEMP 302.3 
LW. TEMP 302. 1 
FL. TEMP 300.2 
EMS(I)= 6. 066E-01 
EMP(I)= 9. 455E-04 
QFC(I)= 8. 225E+00 
0. OOOE+00 
QR(I)= -8. 114E-01 
-1. 405E+00 
Qc(I1I)= -6. 551E+00 
-4. 883E+00 
Pres(kpa) -7. 767E+00 


Upper layer species concentration 


N2 MASS} 2. 61 
PPM ' 1. 894E+05 
co2 MASS !_ 1.3930E-02 
PPM ; 614. 
CO MASS } 3. 678E-04 
PPM 26.7 
H20 MASS : 0.105 
PPM + 1. 190E+04 
OD MASS :  1.316E-04 
1/M : 3. 874E-02 
CT mg-m/1 } 0.155 
URNITURE FIRE SPREADING INFORMATION 
PANEL # IS 1,# OF POLYGONS IS 1 
-RINT FIRE SPREAD VECTORS, MAGNITUDES, AND STATUS 
2. 395982 2. 050000 0. 5000000 0. 1959823 
2. 313682 2. 163682 ©. 5000000 0. 1607716 
2. 200000 2. 246023 0. 5000000 0. 1960234 
2. 100000 2. 148622 0. 5000000 0. 1405170 
2. 100000 1. 952810 ©. 5000000 0. 1397684 
2. 150000 1. 951405 0. 5000000 0. 1198841 
2. 200000 1. 950000 0. 5000000 9. 9999860E-02 
2. 296798 1. 950000 0. 5000000 0. 1396075 
2. 395982 2. 050000 0. 5000000 0. 1959823 
PANEL # IS 2,# OF POLYGONS IS 1 


RINT FIRE SPREAD VECTORS, 


MAGNITUDES, 


AND STATUS 
ee 


ONNOWWOOD 


2. 200000 1. 950000 ©. 5000000 ©. COOOCO0E+00 1 
2. 200000 1. 950000 0. 7000000 0. 2000000 ) 
2. 200000 1. 950000 0. 9000000 ©. 4000000 1 
2. 296798 1. 950000 ©. 9000000 0. 4000000 -1 
2. 296798 1. 950000 ©. 7000000 ©. 2000000 0 
2. 296798 1. 950000 0. 5000000 Q. QOO0000E+00 1 
2. 200000 1. 950000 ©. 5000000 0. COOOC00E+00 1 
PANEL # IS 3,# OF POLYGONS IS 1 
-RINT FIRE SPREAD VECTORS, MAGNITUDES, AND STATUS 
2. 100000 1. 952810 ©. 5000000 0. Q000000E+00 1 
2. 100000 2. 148622 ©. 5000000 0. Q0OO0000E+00 1 
2. 100000 2. 148622 ©. 7000000 0. 2000000 ) 
2. 100000 2. 148622 ©. 7000000 0. 4000000 3 
2. 100000 1. 952810 ©. 7000000 ©. 4000000 3 
2. 100000 1. 952810 0. 7000000 ©. 2000000 0 
2. 100000 1. 952810 ©. 5000000 0. Q0O00000E+00 1 
PRINTING OF ELEMENT STATUS, TEMPERATURE (K), BURNTIME (S), AND BURNRATE (KW) 
1 [eee SI GCOO 45. 17346 3.231312 
2 O 425.6111 ©. QOOOQO0E+00 0. C000000E+00 
3 OQ 430.8394 ©. QOCOO0000E+00 0. C0O00000E+00 
4 OQ 392. 6586 ©. COOQOO0E+00 0. COCO000E+00 
3 1 711.3000 20. 79620 2. 442346 
6 © 420.7428 OQ. OCOO000E+00 0. 0000000E+00 
7 it ico 21. 57066 2. 640063 
=} OQ 444. 7470 ©. OOOOCOCE+00 0. 0O000000E+00 
9 jk eR efolole: 21. 44027 2. 607153 
10 O 415.5579 ©. COOD000E+00 0. 0000000E+00 
11 ey SH elevate: 20. 74713 2. 429640 
io © 393. 6865 Q. OOCOQ000E+00 0. 0000000E+00 
13 OQ 599. 6725 ©. DCOOOC00E+00 0. CO00000E+00 
14 0 567. 4828 ©. DCOOOO00E+00 0. CO00000E+00 
15 0 592.5883 ©. QOOOCO0E+00 0. C000000E+00 
16 OQ 576. 4911 0. COOD000E+00 0. 0000000E+00 
FURNITURE BURNRATE IS 13. 09751 KW AT TIME 70. 00000 SEC 


SMFAC < 0.5 IN RMTERM SUBROUTINE MEANS AO. 3739495 
RELATIVE CHANGE IN THE SUM OF SURFACE-GAS EXCHANGE AREAS 
PROVIDED THE SUM OF GAS-GAS EXCHANGE AREAS ARE CORRECT: 

THE SUMS WERE FOR THE GAS VOLUME # 


1 
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